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Abstract 

We define the edge reconnecting model, a random multigraph evolving in time. 
At each time step we change one endpoint of a uniformly chosen edge: the new 
endpoint is chosen by linear preferential attachment. We consider a sequence of 
edge reconnecting models where the sequence of initial multigraphs is convergent 
in a sense which is a natural generalization of the notion of convergence of dense 
graph sequences (see [ID])- We investigate how the limit object evolves under the 
edge reconnecting dynamics if we rescale time properly: we give the complete char- 
acterization of the time evolution of the limit object from its initial state up to 
the stationary state using the theory of exchangeable arrays, queuing and diffu- 
sion processes. The number of parallel edges and the degrees evolve on different 
timescales and because of this the model exhibits subaging. 



1 Introduction 

We introduce the edge reconnecting model, a random multigraph (undirected graph with 
multiple and loop edges) evolving in time. Denote the multigraph at time T by Qn{T), 
where T = 0,1,2,... and n = |l^(^„(r))| is the number of vertices. We denote by 
m = \E{Qn{T))\ the number of edges (the number of vertices and edges does not change 
over time). Given the multigraph Qn{T) we get ^„(T+1) by uniformly choosing an edge in 
E{Qn{T)), choosing one of the endpoints of that edge with a coin flip and reconnecting the 
edge to a new endpoint which is chosen using the rule of linear preferential attachment: 
a vertex v is chosen with probability 2m+~nK ' where d{v) is the degree of vertex v in QniT) 
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and K G (0, +00) is a fixed parameter of tlie model. A similar dynamical random graph 
model is studied in 

Our aim is to describe the time evolution of the edge reconnecting model Qn{T) 
when 1 ^ n using the terminology of dense graph limits. The notion of convergence of 
simple graph sequences was defined and several equivalent characterizations of graphons 
(limit objects of convergent simple graph sequences) were given in [TU]. In [H] we give a 
natural generalization of the theory of dense graph limits to multigraphs (see also [TT] 
for similar results in a more general setting). Denote by A4 the set of multigraphs. We 
say that a sequence of multigraphs {Gn)'^=i is convergent if for every F G the limit 
g{F) = limn^oo 't={F, Gn) exists, where 

t^{F,G) = -^ J2 n'^v,weV{F): E{v,w) = Ei^{v),ipiw))] (1) 

and E{v, w) denotes the number of edges between v and w in F (loop edges count twice), 
moreover g{-) is a "non-defective probability distribution" on the set of multigraphs (see 
Subsection 12.21 for details). t={F,G) may be interpreted as the density of copies of F in 
G, thus the sequence {Gn)^_^ is convergent if the density of every fixed graph F in Gn 
converges as n — )■ cxd, and "no mass escapes to infinity" during this limiting procedure. 

The definition of the limit objects of convergent multigraph sequences is slightly more 
complicated than that of graphons. A measurable function W : [0, 1] x [0, 1] x No — > [0, 1] 
satisfying 

00 

W{x,y,k) = W{y,x,k), ^W{x,y,k) = 1, W{x,x,2k + 1) = (2) 

fc=0 

is called a multigmphon. Note that {W{x, y, k))'^^^ is a probability distribution on No for 
each x,y E [0, 1]. We say that Gn ^ W if for every F G we have lim„^oo t={F, Gn) = 
t={F,W) where 

t={F,W):= T\ W{x^,x^,E{v,w))dxidx2 ... dxk. (3) 

Theorem 1. of [Ij states that if a sequence of multigraphs (G„)^^ is convergent then 
Gn — W for some multigraphon W and conversely, every multigraphon W arises this 
way. We say that a sequence of random multigraphs (Gn)^_^ converges in probability 

to a multigraphon W (or briefly write Qn W) if for every simple graph F we have 
t={F,gn)^t={F,W),i.e. 

VFGA^Ve>0: lim P (|t=(F, ^„) - t=(F, W)| > e) = 0. (4) 

In Lemma 2.1 of pjij we build on methods and models of Section 3.4 of [1] to explicitly 
describe the unique stationary distribution ^„(oo) of the edge reconnecting model and 
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in Theorem 2 of [12] we prove that there is a multigraphon Woo such that 

Qn{00) Woo, n ^ OO (5) 

under the condition that m ^ ^pn^, where p G (0, +oo) is a fixed parameter of the model 
called the edge density. The form of the limiting multigraphon Woo depends on p and 
the linear preferential attachment parameter k. 

Now we describe the main results of this paper: if we consider a sequence of edge 
reconnecting models with a convergent sequence of initial multigraphs ^„(0) — )■ W (sat- 
isfying some extra regularity conditions), then for every t G (0, +oo) we have 

Qnit ■ n^) A W and G^t ■ n^) ^ W, n ^ oo (6) 

where the multigraphons Wt and Wt are explicit, continuous functions of t, the initial 
multigraphon W and k. Moreover we have 

lim Wt = W, hm W = lim Wt, lim Wt = Woo (7) 

where Woo is the multigraphon in Thus by ([7]) the convergence theorems 1^ give 
the full characterization of the time evolution of the multigraphons arising as the graph 
limits of the edge reconnecting model. 

Although our theorems are stated using the "multigraphon" formalism, in their proofs 
we use the correspondence between the theory of graph limits and that of exchangeable 
arrays, a connection first observed in [B] . The basic idea of the proof of our main theorems 
is to relate the time evolution of the edge reconnecting model to certain continuous-time 
stochastic processes using an appropriate rescaling of time: 

• If we fix a vertex v G V{Qn{0)) and denote by d(T,v) the degree of v in Qn(T) 
then the evolution of the R+-valued continuous-time stochastic process jid{n^ ■ t, v) 
"almost looks like" that of a Cox-Ingersoll-Ross process (a diffusion process that is 
commonly used in financial mathemathics to model the evolution of interest rates). 
This fact is rigorously proved using the theory of stochastic differential equations 
and is used in the proof of Qnit ■ 'n?) Wt. 

• If we fix two vertices G V^(^n(0)) and denote by E(T,v,w) the number of 
parallel/loop edges connecting v and w in Qn(T) then the evolution of the No- 
valued continuous-time stochastic process E{n^ ■ t,v,w) "almost looks like" that 
of the queue length of an M/M/oo-queue. This fact is rigorously proved using a 

coupling argument and is used in the proof of Qn{t ■ n^) Wt- 

The most interesting property of the edge reconnecting model is the separation of two 
different timescales in ^ and ([7]): the degrees of the vertices only change significantly 
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on the re" timescale, whereas the number of parallel (or loop) edges between two vertices 
evolves on the much faster 'n? timescale. The arrival rate of the M/M/oo-queue describing 
the evolution of E{n'^ ■ t, v, w) depends on the current degrees of v and w (if their degrees 
are high then edges appear between them with higher probability due to preferential 
attachment), but since the degrees evolve on the much slower timescale, they may be 
treated as constant background parameters on the n"^ timescale. The stochastic process 
E{n^ ■ t + ■ s,v,w) looks stationary in the time variable s G M if t G (0, +00) is fixed 
and 1 <^ n, but different values of t yield distinct pseudo-stationary distributions since 
■ {h — ii) steps are enough for the background variables (degrees) to significantly 
change. This phenomenon is called subaging in [Ij. 

The rest of this paper is organized as follows: 

In Section |2] we introduce some notation, precisely formulate the above stated results, 
and give some heuristic hints on their proofs. 

In Section |3] we relate the theory of multigraph limits to exchangeable arrays. 

In Section H] we prove some technical lemmas showing that degrees and multiple edges 
in the edge reconnecting model are well-behaved. 

In Section |5] we prove the rigorous version of Gnit ■ n^) Wt, Theorem m 

In Section |6] we prove the rigorous version of Gnit ■ n^) ^ Wt, Theorem 12 

Acknowledgement. The author thanks Laszlo Lovasz and Yannis Kevrekidis for posing 
the research problem that became the subject of this paper. The comments of the 
anonymous referees helped a lot in the developement of the paper to its current form. 

The research of Balazs Rath was partially supported by the OTKA (Hungarian Na- 
tional Research Fund) grants K 60708 and CNK 77778, Morgan Stanley Analytics Bu- 
dapest and Collegium Budapest. 

2 Notations, definitions, theorems 

This section is organized as follows: 

In Subsection 12.11 we precisely define the edge reconnecting model. 

In Subsection 12.21 we give a probabilistic meaning to t={F,W) by introducing W- 
random multigraphs and also define the average degree D{W, x) of W at point x. 

In Subsection 12.31 we recall some relevant properties of the M/M/oo-queue and the 
Cox-Ingersoll-Ross process. 

In Subsection 12.41 we state Theorem [T] and Theorem |51 and we also give some heuristic 
comments on ideas behind their proofs. 

In Subsection 12. 51 we derive ([7]) and relate some properties of the multigraphons from 
our main theorems to the configuration model. 

Denote by No = {0, 1, 2, ... } and [n] :={!,..., n}. Denote by Ai the set of undirected 
multigraphs (graphs with multiple and loop edges) and by the set of multigraphs on 
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n vertices. Let G e A4„. The adjacency matrix of a labeling of the multigraph G with 

[n] is denoted by {B{i,j))^j^-^j where B{i,j) G Nq is the number of edges connecting the 
vertices labeled by i and j. B{i,j) = B{j,i) since the graph is undirected and B{i,i) 
is two times the number of loop edges at vertex i (thus B{i,i) is an even number). An 
unlabeled multigraph is the equivalence class of labeled multigraphs where two labeled 
graphs are equivalent if one can be obtained by relabehng the other. Thus M. is the set 
of these equivalence classes of labeled multigraphs, which are also called isomorphism 
types. We denote the set of adjacency matrices of multigraphs on n nodes by An, thus 

^„ = {SeN^>^" : B^^B,yie [n] 2\B{i,i)}. 

The degree of the vertex labeled by i in G with adjacency matrix B e An is defined 
by d{B,i) :— Yl]=i thus d{B,i) is the number of edge-endpoints at i (loop edges 

count twice). Let m = ^J2^j=i-^{hj) — i X^ILi ^(-^' ^) denote the number of edges. 
Denote by ^™ the set of adjacency matrices on n vertices with m edges. 

We denote a random element of An by X^j. We may associate a random multigraph 
Qn to X„ by taking the isomorphism class of X„. 

We use the standard notation X„ ~ X[j if X„ and Y„ are identically distributed, i.e. 
\/AeAn: P (X„ = A) = P (X; = A). 

If X„ is a random element of An then 

XW:= (X„(z,j))t=i 

is a random element of Ak- 



2.1 The edge reconnecting model 

Now we describe the dynamics of the edge reconnecting model, which is a discrete time 
Markov chain with state space A^: neither the number of vertices, nor the number of 
edges is changed by the dynamics. {X{T,i,j))'^j^^ — X(T') is the state of our Markov 
chain at time T. 

Given the adjacency matrix X(T) we get X(T + 1) in the following way: let k G 
(0, +oo). We choose a random vertex Void{T) with distribution 

p = i I x(r)) = ^^^^^ (8) 

Then we choose a uniform edge EouiT) = {Void{T), W(T)} going out of Void{T): 

P(mr)^i|x(T),v,,(T))^^^|M(^ (9) 

Note that Soid{T) is uniformly distributed over all edges of the graph at time T and given 
^oid{T), Void{T) is uniformly chosen from the endvertices of £oid{T), moreover 

P (>V(r) = ^ I X(T)) = '^^^^^''^ (10) 
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Given X(T), choose Vnew{T) according to the rules of hnear preferential attachment: 
P{VnUT)=^\MT),VoM{T),W{T)) = ^iMTU) + ^ ^^^^ 

Thus Vnew{T) is conditionally independent from VoidiT) and >V(T) given X(T). 
Let£„e«,(T):={V„e«,(T),W(r)}. 

One step of the Markov chain consists of replacing the edge SoidiT) with £new(T): 
X{T+l,t,j) = X{T,t,j) - nVoidiT) = t,W{T) = j] - 1[Vom{T) = j,W{T) = t] + 

nVneUT) = ^, W(T) = j] + l[Vnen,{T) = J,W{T) = z] (12) 

This Markov chain is easily seen to be irreducible and aperiodic on A^. Note that 
for any k < n the N|f'-valued stochastic process ((i(X(T), z))^^^^, T = 0, 1, . . . is itself a 
Markov chain. 

2.2 Multigraphons and VK-random multigraphs 

In this subsection we give a probabilistic meaning to t={F, W) by introducing VF-random 
multigraphs and also define the average degree D(W, x) of W at point x. Note that the 
notion of the VF-random graph (see Definition 12. ip is already present in [lOj . 

Suppose F E Aik, G E A4n and denote hj A E Ak and B G An the adjacency 
matrices of F and G. If g : Ai ^ M. then we say that is a multigraph parameter. Let 
g{A) := g{F). Conversely, if g : IJ^^^ — t- M is constant on isomorphism classes, then 
g defines a multigraph parameter. 

We define the induced homomorphism density of F into G by 

UF,G):=UA,B):=^^ J2 l[yi,je[k]:A{z,j) = B{^{2),^{j))]. (13) 

ip:[k]^[n] 

We say that a sequence of multigraphs (Gn)^^ is convergent if for every k E N and 
every multigraph F E Aik the limit g{F) = limn^oo t={F, Gn) exists, moreover we have 
J^AeAk ~ ^'^^ every multigraphon W (see ([2])) and multigraph F E Aik with 
adjacency matrix ^4 G we define 



t={F,W) := t={A,W) := 




We say that Gn ^ W if for every F G we have limn^oo t={F,Gn) = t={F,W). 
By Theorem 1. of [9] we have that a sequence of multigraphs {Gn)^^i is convergent then 
Gn W for some multigraphon W. The limiting multigraphon of a convergent sequence 
is not unique, but if we define the equivalence relation 

Wi = W2 ^ VFeM: t={F, Wi) = t={F, W2) (15) 



6 



then obviously Gn Wi <^=^ G„ — )■ W2. For other characterisations of the equivalence 
relation = for graphons, see [3J. 

If X„ is a random element of An for each n E N, Qn is the isomorphism class of X„ 
and is a multigraphon, then we say that X„ W if Qn W, see @. 

For a multigraphon W and a; G [0, 1] we define the average degree of at x and the 
edge density of W by 

,1 00 



D{W,x):= / ^A;-iy(x,?/,fc)d|/ (16) 

„i „i 00 

p{W):= / / ^fc- H^(x,|/,A;)d?/dx (17) 
Jo Jo f^^o 

If p{W) < +00 then D(VF, x) < +00 for Lebesgue-almost all x. 

We say that a [0, 1] -valued random variable U is uniformly distributed on [0, 1] (or 
briefly denote U ~ U[0, 1]) if P ([/ < x) = x for all x G [0, 1]. 

Definition 2.1 (ly-random multigraphons) . 

Fix k eN. Let {Uifl^^ he i.i.d., Ui ~ [/[0, 1]. Given a multigraphon W we define the 
Ak-valued random variable = {Xv/{i, .^^ as follows: 

Given the background variables {Ui)\^^ the random variables {X]y(i, j)).^^^^^ are con- 
ditionally independent and P (^Xw(i,j) = I [ ([/j)jL^) = W{Ui, Uj, I), that is 

WAeAk-. P (xiJ) = A I ([/,)ti) := llW{U„U„A{z,j)). (18) 

i<j<k 

In plain words: if i 7^ j and Ui = x, Uj = y then the number of multiple edges 
between the vertices labeled by i and j in X^y has distribution (W{x,y,l)^'^^ and the 
number of loop edges at vertex i has distribution (W{x, x, 2/))^^^. 

For every multigraphon W and A G Ak we have 

t4Aw^)™=°^p(xl5 = A). (19) 

Recalling it follows that Wi ^ W2 if and only if VA; G N : xJJ], ~ xjj^, thus 
the distribution of the VF-random multigraphons determine the multigraphon up to = 
equivalence. Recalling f|T6l) and f|T7|l we have 



D{W, x) = E {Xw{l, 2)\Ui = x), p{W) = E (X^(l, 2)) . (20) 

Note that the weak law of large numbers (heuristically) implies that 

^-d{X^;i,i)^D{W,U,), l«n (21) 



n 



This relation is the reason why we gave the name average degree to D{W, x) 
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2.3 Auxiliary stochastic processes 

In this subsection we recall the definition and some properties of two stochastic processes: 
the M/M/oo-queue and the C.I.R. process. 

First recall the formulas defining the Poisson, binomial and gamma distributions: 

p(fc,A):=e-^^ (22) 



h{k,n,p):= \^^jp^{l-p)--^ (23) 

Oct — 

g{x,a,(3) := l[x > 0] (24) 

We say that a nonnegative integer-valued random variable X has Poisson distribution 
with parameter A (or briefiy denote X ~ POI (A)) if P {X = k) = p{k, A) for all k E N. 
We say that a [n]-valued random variable Y has binomial distribution with parameters 
n and p (or briefiy denote Y ~ BIN(n,p)) if P [Y = k) = h{k,n,p) for all k G [n]. 
We say that a nonnegative real-valued random variable Z has gamma distribution with 
parameters a and (3 (or briefiy denote Z ~ Gamma(a, /3)) ifP{Z<z) = g(x, a, /3) dx. 

The M/M/cxD-queue with arrival rate n and service rate 1 is an Nq- valued continuous- 
time Markov chain Yt, t e [0, +oo) with infinitesimal jump rates 

(25) 
(26) 
(27) 

Heuristically, Yt is the length of a queue at time f, where customers arrive according to a 
Poisson process with rate /i, customers are served parallelly and each customer is served 
with rate 1. It is well-known (see Exercise 5.8 of [5]) that if Iq = ^ ^ ^^o then 

k 

P{Yt = k\Y^ = h)= q(t, h,k,^l)■=Y, b(/, K e-*) ■ v{k - /, (1 - e-*)/i) , (28) 

z=o 

i.e. Yt has the same distribution as the sum of two independent random variables with 
BIN(/i, e~*) and P0I((1 — e~*)yu) distributions. From fl28|) we get that indeed Yt k 
as t — )■ and the stationary distribution of the queue is POI (/i): 

limq(t, /i, A;, /i) = IL[A; = /;,], lim q(t, /i, fc, yu) = p(A;, /i). (29) 

Fix G (0,+oo). The Cox— Ingersoll— Ross (C.I.R.) process is a diffusion process 
with stochastic differential equation 



P [Yt+At = k + l 




= k) 


= ljdt + o{dt) 


P {Yt+^t = k-l 




= k) 


= kdt + o{dt) 


P {Yt+At = k 




= k) 


= 1 - (/i + k)dt + o{dt) 



dZt ={^^- dt + \/2ZtdBt 



(30) 
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where Bt denotes the standard Brownian motion (for an introduction to SDE, see 

Heuristically the SDE f l30p tells us the mean and variance of small incerements of the 
continuous-time ]R+-valued Markov process {Zt)^^^ given the present value of Zf. 

E [Zt+at -z\Zt = z) ^ (^K-'^z^ dt, [Zt+dt - z\Zt = z) ^2zdt (31) 

It is well-known (see Chapter 4.6 of ^) that if we denote 
a := — and r(a, t) :- 



p ' exp(at) — 1 

and if we start the process (Zj)^>q from the initial value Zo = z then 2{T{a,t) + a) ■ Zt 
follows a noncentral chi-square distribution with 2k degrees of freedom and non-centrality 
parameter 2z ■ r(a, t), thus we have P (Z^ < x \ Zq = z^ = f(t, z, y) dy where 

oo 

/(^: z^y) = ^ P(«, z ■ r(a, t))g{y, n + i, r(a, t) + a). (32) 

Note that using fl32l) one can derive that indeed Zt — ^ 2; as t — t- and the stationary 
distribution of (Zj)^>q is Gamma(K, ^): 

\im fit, z,y) = 5:,^y, Urn f{t,z,y) = g{y,K,-). (33) 
2.4 Statements of Theorem [1] and Theorem [2] 

In this subsection we state the main results of this paper describing the time evolution 
of the limiting multigraphons of a sequence of edge reconnecting models X„(-), n — > 00. 

In Theorem [1] we precisely formulate X„(t ■ n^) — ^ Wt- 

In Theorem [2] we precisely formulate X„(t ■ n^) — ^ Wt. 

Note that in ([6]), ([7j) and above we used the notations Wt and Wt in order to give the 
most simple formulations of these results, nevertheless our real notations are going to be 
slightly different. 

Now we describe the evolution of the edge reconnecting model by describing the 
evolution of the limiting multigraphons. We consider a sequence of initial multigraphs 
(^n)^i which converge to a multigraphon W. We assume |l^(G'„)| = n. We denote the 
adjacency matrix of by i?„ G An- We assume that the technical condition 

1 " 

3A>0, C<+cx)Vn: Yl e^^-^^'^'^ < C, - ^ e^^-^^-^^ < C (34) 

\2/ i<j<n 1=1 

holds. 
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First we state Theorem [T] about the evolution of the edge reconnecting model on the 

T = 0{n'^) timescale. In the Introduction this result was referred to as Qn{t ■ n^) Wt 
in order to make the notation as simple as possible. In fact we are going to prove 
Gn{t ■ ■ n^) Wt, thus Wt = Wtf The notation Wt will no longer be used. 

Theorem 1. Let us fix k E (0,+oo). We consider the edge reconnecting model ^n{T), 
T = 0, 1, . . . on the state space An^"^^ and initial state X„(0) = G A^^"^^ for n = 
1,2,.... We assume Bn ^ W for some multigraphon W and that flM|) holds. 
Then for all t G [0, +oo) we have 




^ p{W) ■ r? 



Wt as n oo (35) 



1[2 I k] ■ EZo y, hMt, I, I, ^fx = y ^ > 

We prove Theorem [T] in Section O Before stating further theorems, we devote a few 
paragraphs to the heuristics behind Theorem [H 

In order to give some insight about (1361) . we now give a probabilistic way to generate 



a random element of Ak with the same distribution as x|^^ (see Definition 12. ip : 

We first generate using the background variables (f/j)f=i, then we get x|^^ by 
letting the entries (^Wtih j=i evolve in time: 

• if z < j, we run an M/M/oo-queue Yt with initial value Yq = ^w{hj), arrival rate 
^^^'^p{W)^'^^^ ^^"-^ service rate 1 and let 'Kw^{i,j) := Yt 

• if z = J, we do the same thing with the only exception being that the queue 

D{W,UiyD{W,Uj) 



describing the evolution of the number of loop edges has arrival rate 



2p(W) 



Now we give a heuristic argument explaining why do M/M/oo-queues enter the pic- 
ture: 

We look at the evolution of Xn(T,i,j) for some i ^ j (the case of loop edges is 
analogous). We denote by VniT,i) := id(X„(T),i). From (pTO]) and it follows that 



P (X„(T + 1,.,,) = X„(r,.,,) + 1 I X„) ^ 2m + nK + 

Vn{T,j)-n Vn{T,i)-n + K 1 P„(r, 2)D„(r, j) 



2m 2m + hk m p 



(37) 



P (X„(T + l,z, j) = Xn{T,i,3) - 1 I X„) ^ -X^{T,i,j) (38) 
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n2 



In the statement of Theorem [T] we used the time scahng T = [t ■ ^ t ■ m, thus if 
we denote dt := — then T + 1 corresponds to t + dt. If we define Yt := X{t ■ m, i, j) and 
compare fl37|) . fl38|l to fl25l) . fl26l) then we see that the time evolution of Yt approximates 
that of an M/M/oo-queue with arrival rate /i = ^n^j service rate 1. We will 

later see that on the T x timescale T>n{T, i) does not change significantly, so that we 
have 

We note here that the identity D{W, x) = D{Wt, x) can be formally derived from (15^ . 

Now we look at the evolution of the edge reconnecting model on the T = 0{n^) 

timescale. In the Introduction this result was referred to as Qn{t ■ n^) Wt in order to 

make the notation as simple as possible. In fact we are going to prove Qn{t- p{W)-n^) — ^ 
Wt, thus Wt = Wpt- The notation Wt will no longer be used. 

Theorem 2. Let us fix k, & (0, +oo). We consider the edge reconnecting model X„(T), 
T = 0, 1, . . . on the state space An^^^ and initial state X„(0) = Bn E ''"^ for n = 
1,2,.... We assume ^W for some multigraphon W and that fl34|) holds. 
Then for all t G (0, +oo) (but not for t=0) we have 

^n{lt- p{W) -n^l) ^Wt as n^oo (39) 

where 

and F^^ is the inverse function of Ft{x) = f(t, y) dy where (recall (132|) and note that 
f{t,z,y) also depends on the parameters n and p) 

POO 

f{t,y)= / f{t,z,y)dF,{z), (41) 



and Fq{x) = l[D{W,y) < x] dy, x G [0, +oo). 

We prove Theorem [2] in Section [HI Now we devote a few paragraphs to the heuristics 
behind Theorem [T] In order to give some insight about f l40p . we now give a probabilistic 
way to generate a random element of Ak with the same distribution as : 

Let us first generate using the background variables (Ui)'^^-^. Let us define 
Zi{0) := D(W,Ui). Now (-^^(0))^^^^ are i.i.d. with probability distribution function Fq. 
We let (Zj(t))^>Q evolve in time according to the SDE fl5Ul) . so that (^i(t))f=i are i.i.d. 
with probability distribution function Ft. Given the background variables {Zi{t))^^-^ let 
X^^(z, j) ~ POI(^^^) if z ^ J, and let iX^^(z, t) ~ POI(^0|ff ). 
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Now we give a heuristic argument explaining why do C.I.R. processes enter the picture: 
Pick i e [n] and denote by V^iT) := id(X„(T),z). It follows from dS} and that 

E {V4T + 1) - I?„(T) I X„(T)) = ^ " - ^-^(^- -V„iT)] (42) 

^ ' 2m + UK 2m 2mn \ p J 

D= (i>„(r + 1) - v.m I x„(r)) . ^^^^ + ^ . — 2i>„(r) (43) 

We look at the time evolution of the stochastic process Zf := Vn{\t ■ 2nm\). In the 
statement of Theorem |2] we used the time scaling T = [t ■ p ■ ~ t ■ 2mn. If we let 
dt = then T + 1 corresponds to t + dt. Let dZt := Zt+dt — Zt- From ( l42l) and (H3l) 
we get 

E (dZi I Zt) ^ (^«: - -^Z^ dt (dZt I Zt) ^ 2Ztdt 

Thus the process Zt approximates the solution of the SDK of the C.I.R. process (13T]) . 

2.5 Properties of Wt and 

In this subsection we relate some properties of the multigraphons Wt and Wt that appear 
in Theorem [1] and Theorem [2] to the configuration model. But first, we show 

lim Wt = W, lim Wt = lim W, lim W = W^ (44) 

t— >0+ t— >oo t— >0+ i— >oo 

where Woo is the multigraphon defined by 

F'H^)P-Hy) ^ if a; ^ y 



and F ^ is the inverse function of F{x) = g{y, k, ^)dy, see 

In order to make sense of ( H4|) we define convergence on the space of multigraphons: 
we say that lim„^oo Wn = W ii lim„^oo t={F, Wn) t={F, W) for all F e A^. The hmit 
of a convergent sequence is only determined up to the = equivalence, see ( fT5|) . 

One can see by looking at f l2^ . (12^ that IVt is a continuous function of t. 

Similarly, (gn]), dH]), (132]), (122]), (jM]) imply that Wt is a continuous function of t. 

If we substitute t 0+ into (|36]) we indeed get W7 W bv fl29|) . 

If we substitute t oo into (|40]) we get W Woo by fl33|l . 

If we let t ^ oo in ([36]) and t ^ 0+ in (HQ]) we get 



hm Wt(x,2/,A:)^^ A VW^T^M^ -f ^^^^ 
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where Fq(z) = '5\.[D(W,y) < z]dy and Fq ^(x) := min{2; : Fq(z) > x}. We have 
hmt_j.oo Wt = hmt^o+ Wt, because the corresponding ly-random muhigraphons have the 
same distribution, since D{W, Ui) ~ FQ^{Ui). Thus we have seen that (jH]) holds. 

A well-known way to generate a random multigraph with a prescribed degree sequence 
is called the configuration model: we draw d{v) stubs (half-edges) at each vertex v and 
then we uniformly choose one from the the set of possible matchings of these stubs. 
In [13] we call such random multigraphs edge stationary and in Theorem 1 of [13] we 
characterize the special form of limiting multigraphons that arise as the limit of edge 
stationary dense multigraph sequences: these multigraphons are of form fH5l) where F is 
a generic probability distribution function on M_|_ and F~^{x) := min{2; : F{z) > x} is 
the generalized inverse of F. The name of edge stationarity comes from the fact that the 
space of edge stationary distributions is invariant under the edge reconnecting dynamics, 
see Section 4 of p^. The stationary distribution of the edge reconnecting model is an 
example of an edge stationary multigraph, see (l5l). fH5l) . 

The heuristic explanation of the fact that the limiting multigraphon lim(_j.oo Wt from 
(H^ has the special form that appears in Theorem 1 of [13j is as follows: If T ^ t ■ 
where 1 -C t then the degrees of vertices in ^n(O) and Gn{T) are very close to each other 
(c.f. D{W,x) = D{Wt,x)), whereas T steps are enough for the model to rearrange and 
mix the edges, so Gn{T) looks edge stationary. Similarly, for any t > 0, Wt from fj40l) 
also looks edge stationary. 

Roughly speaking, if we start the edge reconnecting model from an arbitrary initial 
multigraph, then we have to run our process for <^T steps until QniT) becomes "edge 
stationary" and run it for n? <^T steps until QniT) becomes "stationary". 

3 Vertex exchangeable random adjacency matrices 

In this section we define the notion of vertex exchangeability of random adjacency matri- 
ces and recall two lemmas from pTS]: in Lemma [3?T] we relate convergence of dense random 
multigraphs to convergence of the probability measures of the corresponding vertex ex- 
changeable random arrays and in Lemma 13.21 we give sufficient conditions under which 
convergence of dense random multigraphs imply convergence of the degree distribution 
of these graphs, see fl2Tl) . 

Let X = (^(«, denote a random element of An- We say that the distribution 

X is vertex exchangeable if for all permutations r : [ra] — [n] and B e An 

P(V2,jgM: B{i,3) = X{i,3))='P(yi,3e[n]: B{i, j) = X{t{^,t{j))) , 

that is {X{i, j))^ and {X{r{i),T{j)))'^ have the same distribution: 
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In graph theoretic terms fl48l) means that the distribution of the random graph is invariant 
under the relabehng. It follows from Definition 12.11 that is vertex exchangeable. 

In the statements of Theorem [1] and Theorem [2] the initial state of the Markov chain 
X„(T) = {Xn(T, was the deterministic adjacency matrix X„(0) = but if we 

define 

X„(0,^,J):=5„(7r(^),7r(J)). (49) 

where tt denotes a uniformly chosen random permutation of [n] and denote the edge 
reconnecting Markov chain with this initial distribution by X„(T), T = 1,2,..., then 

(Xn{T,t,j)y ~ (X„(r,7r(^),7r0-)))" , t=(A X„) ~ t=(A Xn), (50) 

thus we get that the assertion of Theorem [1] and Theorem [2] holds for X„(T) if and only 
if it holds for X„(r). From now on we are going to use this trick to replace X„(T) by 
X„(T) and assume that the distribution of X„(r) is vertex exchangeable. 

If X„ is a random element of An for each tt, G N and is a multigraphon then we 
say that X„, — ^ X^y if for all A; G N we have xL^' x{^, i.e. 

V A; G N, A G A : hm P (x|fl = A) = P (x\S = a] ^ t={A, W) (51) 

Recall that we say that X„ if 

Vfc G N VA G A V£ > : lim P(|t=(A,X„) -t=(v4,W^)| > e) = 0. (52) 

n—>OD 

We state here Lemma 3.1 of [13j without proof: 

Lemma 3.1. Let X„ = {Xn{i, be a random, vertex exchangeable element of An 

for all n ^N. The following statements are equivalent: 

Xn^W ^ Xn^Xw (53) 

For a real- valued nonnegative random variable X define E (X; m) := E (X ■ l.[X > m]) . 
A sequence of real-valued nonnegative random variables {Xn)'^^^ is uniformly integrable 
(see Chapter 13 of [H]) if 

lim maxE (X„; m) = 0. (54) 

m—^OD n 

We state here a special case of Lemma 3.2/(ii) of without proof: 

Lemma 3.2. // X„ is a random vertex exchangeable element of An for each n G N, 

X„ Xvi/ holds for some multigraphon W and the sequences (X„(l, 1))^]^ and (X„(l, 2)) 
are uniformly integrable then for all k we have 




where Xy- is generated using the background variables {Ui)^^^ according to Definition \2.1\ 
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4 Bounds on multiple edges and degrees 

In this section we state and prove Lemma l4m which, roughly speaking, states that degrees 
and multiple edges in the edge reconnecting model remain well-behaved. 

If we replace the initial matrix Bn with its vertex exchangeable version (^^(O, i, 
using the trick fHU]) then the technical condition f l5^ becomes 

3 A > 0, C < +00 Vn Vz,i G [n] : E (e^^"(o«)) < C (56) 

It is easy to see that fISB]) implies that the sequences (X„(0, 1, 2))^^^ and (X„(0, 1, 1))^^ 
are uniformly integrable. If we assume Bn — > W and define p := p{W) then 

lim ^ = hm 1 X: ^r.iO,^,J) = lim 1 ±EiX40,^,J)) ^ 

i,j=l i,j=l 

hm f^E(X„(0,l,2)) + lE(X„(0,l,l))') E (Xh.(1, 2)) # p (57) 

Now we state and prove a lemma which says that if the initial state X„(0) of the edge 
reconnecting model is well-behaved (i.e. (156|) holds) then model remains well-behaved at 
later times T as well: 

(i) For all T = 0{n^) the normalized degree D(T) = ic/(X„(T),i) of a vertex i G [n] 
satisfies D(T) = (9(1) uniformly in n, a bit more precisely: P (-D(T) > z) decays 
exponentially as ^ — )■ 00. 

(ii) For all T = 0{n^) the number of parallel/loop edges Xn(T,i,j) between vertices 
i,j G [n] satisfies X„(T, = 0{1) uniformly in n, a bit more precisely: X„(T, 
has finite moments. 

(iii) if Ti < T2 = 0{n^) and T2 — Ti ^ then D(Ti) ^ D(T2), a bit more precisely: 
the second moment of D(Ti) — D(T2) is 0{n^^{T2 — Ti)). 

Lemma 4.1. Let us fix k, p & (0, +00). We consider the edge reconnecting model X„(T), 
T = 0, 1, . . . on the state space A^^"'^ with a vertex exchangeable initial state X„(0) for 
n = 1,2,... satisfying lim„_>.oo ^^2"^* = P o-nd (13^ for some A < ^. Then 

(i) There exists an n' & N such that for every t, z & [0, +00) and n > n' we have 

P ( max -(i(X„(r),i) > z] <C- e^^^' ■ e"^" (58) 

\0<T<2mntn ^ ^' ^ - ^ - ^ ^ 

With the C 0/ dSHD. 
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(ii) For every p > 1 and t G [0, +00) there exists a C = C'{k, p, A, C,p, t) (where C is 
the constant from fl56p ) such that for all n E 'N, i,j E [n] and T < 2mnt we have 

E(X„(T,^,jr)<C'- (59) 

(Hi) There exists a constant C" = C"{k, p, A, C, t) such that for all n eN, all Ti < T2 < 
t ■ and i G [n] we have 

E((-d{X4T2),t)--d{X4Tr),2)] I < C"-{T2-T,) 



Proof of Lemma \4-l\ (jl]) • Fix i G [n] and denote 



diT) := d(X„(^),^), D{T) := -d(X„(T),z). 

n 

Denote by (^r)o<T natural filtration generated by the process. 

If a{T) := E"(e^-(^(^+i)-^(^)) - 1 | ^t) then M{T) := e^^(^) nr=o'(l + ^(O)"' is a 
nonnegative martingale. By the submartingale inequality we have 

P max M(T) >x]< ^ (61) 

\0<T<T' J X 

^max ^ M(r) <x =^ VT < T' : e^^^^^ < xexp j^J^ a(/)^ (62) 
Now we give an upper bound on a{T). Using 



, 1=0 



D{T + 1) = D{T) + ^llVne^T) =t]- ^l[VoUT) = t], (63) 

dH]), dUD and the fact that Vnew(T) and VoidiT) are conditionally independent given J^t 
we get 

= f 1 + ^(I)±^(ei - 1)) f 1 + 'J^(e-^. -,))-!< 
(i(r) + « /^A _^ ie-— ^ + (-- + -e- 



2m + nK \n 2 n^y 2m \ n 2 ra' 
^ ^ d{T)- (e^X- (— + Ik) -UK] +2mK- [l + le^-]] (64) 



n + n2mn \ \ \ n 2 J J \ 2 n 
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Now A < - and lim„_>.oo "^2" = P, thus if n is big enough then A < e « ■ — n-r, which 
imphes that the coefficient of d{T) is negative in the right hand side of fl64|) . thus 



1 + i^e^ 

a(T) < ^— A/c < ^— 2Ak. (65) 

^ ^ - 2mn 1 + ^ - 2mn ^ ^ 

2m 



From ([HI]), ([62]) and ([65D it follows that 
P 

Substituting x = exp(— 2KAt) exp(A2;) and using 



/ \ E (e^^^°'>) 
max e^-°(^) > xexp(2AKt) < ^ 

\0<T<2mnt J ^ 



E (e-<«)) = E (^exp (^i ^'^ (». j)) ) £ ^ 



-> exp(AX(0,z,j)) 1 < C 



we arrive at fl58l). □ 



Proof of Lemma \4.1\ ([n]) . Fix n and 2, j G [n]. We only prove the statement of the lemma 
if z 7^ j, the proof of the diagonal case is similar. Denote by 

X{T) := X^{T,z,j), d{T,i) := d(X„(T),z), D{T,i) = -d{T,t). 

n 



Using ([T2]l we get 

P(X(T+1) = X(T) + 1| Jr) 

d{T,i) + K [d{T,j) [ ^ X{T)\\ d{T,j) + K fd{T,i) X{T) 



2m + riK 



f d{T,j) f X{T) \\ d{T,j) + K f d{T,z) f X{T)\\ 
V 2m V d{TJ)JJ 2m + nK V 2m V d{T,t) J J ^ ' 



2m \ 2m + nn J 2m \ 2m + nn 

(67) 

From this it is straightforward to derive 

E (e^^(^+i) - e"^(^) I J^t) < 
gAX(T) j^/ A _ diT,i) + KdiT,j) ^ diT,j) + KdiT,i) \ ^ _ X(r) \ 
\ \ 2m + riK 2m 2m + uk 2m J m J 

Define the stopping time 

. / d{T,i) + Kd{T,j) d{T,j) + Kd{T,i) y 

Ty := mm < i : 1 > — 

1^ 2m + nn 2m 2m + uk 2m m 
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and Xy{T) := X{T)l[ry > T]. Now we prove that for all T G N 

E (e^^-(^)) < max |c , exp{yXe^) (^1 + | • (69) 



It is straightforward to check that 



E (e^^«(^+i) - e'^'y^^^ I J-t) < e'^'y^^^ ((e'-l)^ + (e^^ - 1)^^^ . (70) 
^ ' ^ \ m m J 

If we denote E{T) := E (e^'^"^^^) , take the expectation of flTOjl and use Jensen's inequality 
then we get 

E{T + 1) - E{T) <^({e'- l)y + {e'' - . (n) 

m \ A 



We prove (l69l) by induction. For T = we use (156|) . If E{T) > exp(yAe^), then by ([71 
E{T + 1) < E{T) and if E{T) < exp(?/Ae^) then 

E{T + 1) < E{T) + exp(yAe^)-^^— ^ < exp(?/Ae^) f 1 ' ^""^ ~ 



Having established fl69l) we prove fl59|) by showing that 



oo 



E (X(T)P) < 1 + / P (X(T)P > x) dz < +00. 



P (X(r)^ >x)<V [XyiTY >x) + P {X{T) ^ XyiT)) < _^^^-^ + P (Ty >T) < 

exp(Ax^/P) \T<2nmt J 

Now choosing y = x^^"^^ we indeed get P (X(T)p > x) da; < +oo. □ 



Proof of Lemma \4.1\ (Imj) . Fix i G [n]. We use the notation D{T) = -(i(X„(T), z). We 
say that = C'(^n) if there exists a constant c depending only on K,p,X,C and t such 
that ttn < c ■ bn for all n G N. It follows from (1581) that 



\/T<t-n^: E(D(T)) = 0(1) \/T,T <t-n^ : E {D{T)D{T')) = 0{1). (72) 

Using (l63|) . ([8]), (|TT|) and the fact that Vne«)(r) and Void(r) are conditionally independent 
given J^T we get 

' 2m + TiK 2m 2mn + w^k Am^ + 2mnK 
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E((D(T+1)-D(T))2|J-^) = -^( 



1 fnD{T) + K nD{T) nD{T) + KnD{T)\ 



2m + nn 



+ 



2m 



2m + nn 2m 



)■ 

(74) 



We prove (1^ by induction on T2 — T1. 
E ((D(T2 + 1) - D{T^)f) = E ((D(T2) - D{T^)f) + 

2E (E {D{T2 + 1) - D(T2) I J-T,) (D(T2) - D(Ti))) + E ((D(T2 + 1) - ^(ra))') 

E((D(T2)-D(ri))')+0 



m- 



□ 



We state a lemma about the speed of convergence of the M/M/cxD-queue to its sta- 
tionary distribution. 

Lemma 4.2. Let Yt be an No-valued continuous-time Markov chain with infinitesimal 
jump rates 022]) , ([22]) , and initial state h eNq. Then for allt > and I eNq we have 

P {Yt = /) - hm P {Ys = 1) < -{h + fi) (75) 

s— ^00 

Proof of Lemma \4-S\ According to fl29|) Ys POI(/i) as s — )■ 00. Let 

F/^" ~ BIN(/i, e"*), Yr ~ P0I((1 - e"*)//), y~ ~ POI(e-V) 

be mutually independent random variables. 

By ([28D we have F/^'^ + ~ Yt and + ~ POI(/i). 

p (Ft = /) - lim p = /) = |p (r/'" + Yr = - p (^r' + = I < 

p (r/™ + Y^ Y^ + Y^) < p (r/'" ^ 0) + p 7^ 0) = 

1 - (1 - e~*)^ + (1 - exp(-e-V)) < e~* ■ (/i + /i) 



□ 



5 Proof of Theorem [T] 

In this section we prove Theorem [1] by coupling the evolution of multiple edges between 
the vertices l<i<j<kto {'^^^) independent M/M/ 00- queues. 

Given a random element X of Ak we define the modified adjacency matrix X* in the 
following way: let X*{i,j) := X{i,j) \i i ^ j and X*{i,%) := \X[i,%). 



19 



We assume that the distribution of X„(T) is vertex exchangeable (see the paragraph 
after fl50p ). We are going to prove (155]) using Lemma [3.11 we only need to show that for 
all A; G N and t > we have 

Xjl(Lt-^^^^j) AX15. as „-.oo, (76) 

Note that the evolution of ^Xn'(T), ((i(X„(T), is itself a Markov chain under the 

edge reconnecting dynamics. 

We are going to prove (1761) by coupling the ^^-valued discrete-time process 
to an ^jfc- valued continuous-time process Yn\t) which we define now: 

• The initial states are the same: G [k] : Yn{0,i,j) = Xn{0,i,j). 

• Given Yn\o) = x|f^(0), the evolution of Yn{t,i,j) is a continuous-time Markov 
process for each i,j G [k], the entries (Xnitji, j)),-^j^f^ evolve independently and 

Yn{t,i,j) = Yn{t,j,i), thus Yn\t) is a random element of Ak- 

• The process Y*{t,i,j) is an M/M/cxD-queue (see (l25|) . ( 126|) . ( 127|) ) with service rate 
1 and arrival rate 

^ 2m{n)-{l + l[i = j])- ^ ' 

Now we show that for alH > 

YW(t)AxW as n^oo. (78) 
From the assumptions X„(0) W, ( I56p and Lemma [3.21 it follows that 

Y[fl(0), (^^d{X^{0),^)y J ^ (xl5,(D(l^,t/.))ti) . (79) 

Now (HHD easily follows from this, (^Sil, Definition IS] and ([5B]). 

Denote by Vn(T,i) := ^(i(X„(r), i). We are going to construct a coupling (joint 
realization on the same probability space) of the discrete time ^^-valued Markov chains 

Xn'(r) and yL'^' (m^) for T = 0, 1, . . . such that for any z/ < | we have 



hm P fvT < n'^ : XW(r) = Y^ (^)) = 1- 
n^oo y \m{nj J J 



(80) 



Before proving ([80]) we first assume that it holds and deduce Theorem [T] from it: 

Fix t G (0, +oo). If 2 < zy and n is large enough then It ■ m{n) < rf . It is easy to see 
that dZH]), dHO]) and hm^^oo = p{W) together imply (JTS]). 
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Now we start proving fl80p . 

For [k] we define the matrix Eij G Ak by 

E,^{i',f) := l[z = z', J = /] + l[i= f,j = z'] 

Fix n G N. We introduce the events 

E^{T,^,J) := {xW(r + 1) = X[fl(T) ± i?,,} 



T+ 1 



m 



i?x(T,0) := {xW(r + l)=XW(T)} 



T + 1 



m 



It is straightforward to derive from (128|) that there is an absolute constant (7 such that 
if we define 

k 



then 



P(i?+(r,z,j)|-^T)-^ <Erry(r) 



m 



P{EyiT,t,j)\TT) 



m 



P(£;y(T,0) I J-y) -1 + 



< Erry(r) 



< Erry(T) 



(81) 

(82) 
(83) 

(84) 



From the definition of the edge reconnecting model it follows (similarly to fl66|) and fl67j) ) 
that there is a constant C depending only on k and p such that if we define 



Errx(T) := ^ | 1 + ^ X„(T,z,j) + 5^P„(r,z) ) + 



E- 



d(x„(r),z)rf(x„(T),j ) 

2m- (1 + l[i = i]) 



pi, J 



then 



p{Ej,{T,^,J)\J^T)-^ <Errx(T) 



m 



^*(T,^,J) 



P{E^{T,2,j)\J^t) 
P (i?x(T,0) I J-t) - 1 + ^-<^<^^*(^'^'^') + ^-'^- 



< Errx(T) 

< Errx(T) 



(85) 

(86) 
(87) 
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For any joint realization (coupling) of the discrete time processes XL'=1(T) and Y[fl(^), 
T = 0, 1, . . . define E{T) to be the event that the ^fe- valued increment from T to T + 1 
of these two ^fe-valued processes is the same: 

E(T) := i (i?x(T,0)nEy(T,0))U |J |J {E^AT,^, j) n E^{T,t, j)) 

y ':G{ + ,-} i<j<k 

For any coupling the inclusion 

'T+1 



xW(T) = yW (^^) I n EiT) c |xW(r + 1) = 



m 



^9) 



holds. Let 

Err(T) := 2fc2(Errx(T) + Erry(T)). (90) 

Now if we compare ([H2]) to ([in]), (IHSI) to ([HZI) and ([Hl|) to ([HH]), it easily follows that 
there exists a coupling for which 

P {EiT) I Tr) > ]1[XW(T) = Yjf] (^)] ■ (1 - Err(T)) 

Putting this inequality together with flS^ . multiplying both sides by 

1[VT' < T - 1 : X^^](T') = Y|f] 

and taking the expectation of both sides of the inequality we get 

p (yr <T + i: x|fi(r') = y^ (^^^ ^ > 

p(yT' <T: XW(T') = Yjfl (j^^^ " E (Err(T)) . 
Thus in order to prove (IHOjl we only need to show 

n" 

lim VE(Err(T)) = 0. (91) 



n—^oo • 

T=0 



In the remaning part of this section we prove fl9T|l . 
First we show that if T < n'' then 

E(Errx(T)) = C(n-^/2) . (92) 

Since < | < 3, we have n'^ < nm if n is large enough, thus 

E(P„(T,2)2) = 0(1), E(X„(r,2,jr) =0(1) (93) 
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follow from Lemma [4.11 ([!]) and and Lemma [4. II (jn]), respectively. 



E 



m 



d(x„(r),^)rf(x„(T),j) 



2m- (1 + l[i = j]) 



/i 



o -E(|p„(r,z)p„(r,j) -2^n(o,z)P„(o,j)l) 



72^ 



rO(E(|I)„(T,z)-I)„(0,z)|-I)„(r,j))+E(|I)„(T,j)-^^n(0,j)|-^^n(0,z))) 

^0 ( Ve ((i)„(r, z) - i?„(o, Ve (i)„(r, 

v/E((P„(T,j)-^^n(0,j))')VE(P„(0,02)^ ® 



^1 l^(l) = ^K'^') (94) 



The equation marked by (*) follows from the Cauchy-Schwartz inequality. 
Taking the expectation of and using (1^ . we indeed get 



Now we show that if T < n'^ then 

E(Erry(T)) = O [n"^) 



(95) 



The proof of E (^j'^ii)^^ = C^(l) is similar to that of Lemma 14.11 1^ and we 

omit it, E (/U^j) = 0{1) follows from Lemma [4.11 (ji]). Taking the expectation of f lHT]) we 
get iSD. 



Now if we substitute ([95]) and (|92]) into (|90D we get E(Err(T)) = O (n-^/^) from 
which (PT]) follows using < |. 



6 Proof of Theorem 2 



In this section we prove Theorem [2] in two stages: 

In Subsection 16.11 we prove that the joint evolution of the (normed, rescaled) degrees 
of the vertices 1, 2, . . . , A; behave like independent C.I.R. processes if 1 <^ n. Given this 
result we prove (using the results of Section [5]) that after <^ T steps the state of the 
edge reconnecting model is essentially edge stationary in Subsection 16.21 

6.1 Evolution of degrees 

Lemma 6.1. Let us fix k E (0, +oo). We consider the edge reconnecting model X„(T), 
T = 0, 1, . . . on the state space ^JT''"^ with a vertex exchangeable initial state X„(0) for 
n = 1,2, . . . satisfying and (156|) . We assume X„(0) W for some multigraphon W. 
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Then for all t G [0, +00) and k ^ N we have 



{V^{[t-p{W)-n'\,i)), 



(96) 



where are i.i.d. with distribution function Ft{x) = f{t,y) dy where f{t,x) is 

defined by . 

In order to prove this lemma, we are going to apply a special case of Corollary 2.2 of 
[7], which we reformulate to fit our needs and notation: 

Theorem 3. Let /3 : M — )■ M and 7 : M — t- M &e continuous functions. Assume that the 
stochastic differential equation 



dZt = f3{Zt)dt + 7(^t) dBt 



(97) 



has a unique weak solution with Zq = zq for all zq G R. Let Fo{x) be a probability 
distribution function on M. 

Fix A; G N. For each n let (Vn(T, 'i))jg[fc] discrete time M.^ -valued stochas- 
tic process adapted to the filtration (J-'„ 2-)^^^^. Let 



dV^{T,i) ■.= Vn{T+l,i)-Vn{T,i). 



Suppose 



(r'„(0,i))J^^ {Zo^i)i^^ as 00 



(98) 



where (-^0,4)^=1 are i.i.d. with distribution function Fq. Let m : N — ?■ N and suppose that 
for each t* G [0, +00) and each 1 < i, j < k we have 



sup 

te[o,i*] 



[2-m(n)-n-tJ 



J2 E(rfP,(T,^)| J-„,t) 



1 



T=0 



2m{n) ■ n 



l2m{n)-n-t\ 



T=0 



(99) 



sup 

te[o,t*] 



l2m{n)-n-t\ 

J] Gov {dVn{T, i),dVn{TJ) I J-„,t) - 



T=0 



2m{n) ■ n 



\2m{n)-n-t\ 

^[i = 3]-l\Vn{t,i)) 



T=0 



(100) 



[2m(n)-n-t* J 



E E{{dV^{T,t))H[\dV,,{T,z)\> e]\j^n,T) ^ for all e > (101) 



T=0 

as n 00. 
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Then the distributions of the Ml' -valued continuous-time stochastic processes 

[k],t>0 

converge weakly to the distribution of {Zt^i)^^^^ as n ^ oo in the Skorohod space 



i^), where i>o ^'^s i.i.d. solutions of fl97|) . or briefly: 



iV^ ( L2m(n) ■ n ■ tj , z)),e[.],>o ^ (^*.).eW,*>o (102) 



Proof of Lemma \6.1\ 

We are going to use Theorem[3]to prove that for all k we have (11021) where (Zt^j)^gj^j 
are i.i.d. solutions of f l30|) with initial distribution functions P (Zo.i < x) = -Fo(2^); where 
Fq{x) is defined as in Theorem [2j From this the claim of Lemma lOl indeed follows, since 
by ([57!) we have lim„^oo = p(W^), thus 



(I?„(L2m(n) ■r2-tj,2))^ 



G[A;],t>0 



-(P„(L^-p(M^)-n3j,z)),^, 



G[fc],t>0 



0' 



ie[k],t>0 ' 



ViG[fc],i>0 



are 



from which it follows that for each t > the relation fl96|) holds, where {Zt^ 
i.i.d. solutions of f l30|) . and using f l32|) we get that (Zi_j).gj^j are i.i.d. with distribution 



function Ft{x) = f{t, y) dy where f{t, x) is defined by fj4T|) . 
We need to check that ([98]), ([99]), ffTOOD and ffTOTD holds with 



/3(z) 



K 2;, 



From the assumptions X„(0) - 

(2^n(0,^)),,[,] 



14^, ( 156|) and Lemma [3.21 it follows that 

d 



(D(W^,/7.)),,[,], 

thus by fl2Ul) and the definition of Fq in Theorem[2]we get that the probability distribution 
function of D{W, Ui) is Fq and ^ holds. 

Now we chack that (l99l) holds: 



sup 

tG[0,t*] 



[2-m(?i)-n-tJ 

5^ E(dP„(T,2)| 

T=0 
[2m(n)-n-t* J 



[2m(n)-n-tJ 



T=0 
\2m{n)-n-t* \ 



2m{n) ■ n 



T=0 V / 



2m[n) + riK 2m[n) J 2m[n) ■ n \ p 



1 



r=o 



2m{n) ■ n 



O 



n 



K 



P 2m(n) 



^^n(T,z) + 



mm 



ll8t.l[TT1l.|[63t 

< 



(103) 
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By Lemma O di we have E (r'„(T, i)) = C(l), thus E 
imphes 



— > as n — )■ oo which 



We prove fllOOp by treating the cases i = j and ^ 7^ j separately. 
First we prove f llOOp when i = j. Using ([8]), ( ITT]) . ( 163|) and the fact that V„e?i)(r) and 
Void{T) are conditionally independent given J^n,T we get 



sup 

te[o,t*] 



l2m{n)-n-t\ 



J2 D2(dP„(T,z)|j-„,T) 

[2m(n)-nt*J 



T=0 



2m{n) ■ n 



[2m(n)-n-tJ 
T=0 



< 



T=0 



1 / nVn{T, i) f ^ _ nVn{T,i) ^ 



2m 



nVniT, i) + K _ nVniT, i) + k 



[2m{n)-nt* 

E 



2m + UK 
1 



2m + UK 



2m 
1 



2m{n) ■ n 



2Vn{T, i] 



T=0 



O 



n 



2m[n) ■ n \ \m[n) J \^/ \m[n) 



n 



Vn{T,i)] (104) 



By Lemma O © we have E (P„(T, if) = 0{1), thus E ((HMl) ^ as n -> 00 which 
implies fllOOp for i = j. 

Now we prove fllOOp when i j: 



sup 

te[o,i*] 



[2m(n)-n-tJ 

^ Gov {dV„iT,t),dVr,{T,j) I 

T=0 



< 



[2m(n)-n-i* 

E 

T=0 



2m{n) 2m{n) + n/t 2m{n) 2m{n) + n/t 
\ 2m{n) + riK 2m{n) ) \ 2m{n) + n/t 2m{n) 

[2m(n)-n-t*J 



E 

r=o 



2m{n) ■ n \ \m{n) 



O 



n 



{V^{T,i)+V^{T,j)Y\ (105) 



By LemmaOd) we have E (r'„(T, z)^) = 0(1) and E {Vn{T,jf) = 0{1), which implies 
E (f llOSp ) — 7- as n — )■ 00 which in turn implies f llOOp for i ^ j. fllOip is trivial since 

P(|rf^^n(T,z)|<i)=l. 

Having checked that ( l98|) . (199|) . fllOOp and fllOip holds, we can use Theorem [3] to prove 
that we have (11021) where (Zf^j)^gj^j are i.i.d. solutions of (130|) with initial distribution 
functions Fo(3^)? which finishes the proof of Lemma \6.1\ as described in the beginning of 
the proof. □ 
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6.2 Asymptotic edge-stationarity 

Similarly to Section |5] we assume that the distribution of X„(T) is vertex exchangeable. 
We are going to prove fl39l) using Lemma I3.lt we only need to show that for all A; G N 
and t > we have 

xW(Lt-p(iy)-n3j)^xW. (106) 

Proof of Theorem\^ Let (Z^ j) .^j^j be i.i.d. with distribution function Ft{x) = f{t, y) dy 

where f{t,x) is defined by (HTj) . Define the function p{A, {zi)'^^^) for A & Ak and 
Zi e [0,+cx)), i e [k] by 

By (flSj) and (HUl) . in order to prove (11061) we only need to check that for all A E Ak 

hm P (XW ( [t . p{W) .n'\)=A)=E (p{A, {Z,,)l^)) . (108) 
We (somewhat arbitrarily) fix 2 < z/ < |. Let 



2' 

3 



ro":= Lt-p(w^)-r2^J-KJ. 

It easily follows from < | < 3, Lemma WA\ dm]) and Lemma [6.11 that 

(I)„(To", ^))^,j„ ^ (Z,,)^^[,j as n -> oo. (109) 
Now we couple xL^l(To" + T) to Y^f^ ( ^ ) in a similar fashion as in Section |5l 

• The initial state of Y^f' is G [k] : Yn{0,i,j) = Xr,{T^,i,j). 

• Given X„(T,5^), the entries {Xnit^h i))i<ij<ik evolve independently and Yn{t,i,j) = 
yn{t,j,i). 

• Given X„(rj*), the evolution of Y*(t,i,j) is is an M/M/oo-queue with service rate 
1 and arrival rate 

_ d(X„(To"),2)d(X„(ro"),j) _ V^{T-,i)V^{T^,j) 
^ 2m{n) ■ {1 + l[^ = j]) ^JH^ . + = j])' ^ 



Now we show that 



yW i^]^ XS as n^oo. (Ill] 
myn) 
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First note that 

E(nnp(^-('.^). ,.(friti,i) )) '^E(p,A,z,,tj) ,U2) 

Let := ^[(n)- ^^^n^oo tn = +0C foUows from 2 < u. We have 



P(Y|fl(^n)=A)- limP(YW(.) = A) 



Gil 

< 



exp(-y -5^5]; (E(X„(ro",z,j)) + E ™'P'^exp(-t„)0(l). (113) 

i=l j=i 

Thus ffm]) follows from (1TT2|) and ffTT3ll . 

Using the proof of flSOp we can construct a coupling such that we have 

lim P ( VO < T < n'^ : X^J:HT^ + T)= Y'^'] ( ) ) = 1. 

Now ffTOHD follows from this, + [n"] = [t ■ p{W) ■ and ffTTTjl . □ 
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